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ABSTRACT 

Measurements of double-differential charged pion production cross-sections in interactions of 12 GeV/c 
protons on O2 and N2 thin targets are presented in the kinematic range 0.5 GeV/c < p-^ < 8 GeV/c and 
50 mrad < 6*^ < 250 mrad (in the laboratory frame) and are compared with p-C results. For P-N2 (p- 
O2) interactions the analysis is performed using 38576 (7522) reconstructed secondary pious. The analysis 
uses the beam instrumentation and the forward spectrometer of the HARP experiment at CERN PS. The 
measured cross-sections have a direct impact on the precise calculation of atmospheric neutrino fluxes and 
on the improved reliability of extensive air shower simulations by reducing the uncertainties of hadronic 
interaction models in the low energy range. In particular, the present results allow the common hypothesis 
that p-C data can be used to predict the p-N2 and P-O2 pion production cross-sections to be tested. 
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1 Introduction 



The HARP experiment [T] at the CERN PS was designed to measure hadron yields from a large range of nuclear 
targets and for incident particle momenta from 1.5 GeV/c to 15 GeV/c. This corresponds to a proton momentum 
region of great interest for neutrino beams and far from being covered by earlier dedicated hadroproduction 
experiments (U [3] . The main motivations are the measurement of pion yields for a quantitative design of the 
proton driver of a future neutrino factory [4] , a substantial improvement in the calculation of the atmospheric 
neutrino fluxes ^5] and the measurement of particle yields as input for the flux calculation of accelerator neutrino 
experiments 6;, such as K2K lIllBj, MiniBooNE 9J and SciBooNE [10]. 

Measurements of the double-differential cross-section for tt^ production at large angles by protons in the mo- 
mentum range of 3 GeV/c-12.9 GeV/c impinging on different thin 5% nuclear interaction length (Ai) targets 
have been reported in [TTl [T^l [131 [H] . These measurements are of special interest for target materials used in 
conventional accelerator neutrino beams and in neutrino factory designs. 

The results on the forward production of 7r+ in p-Al interactions at 12.9 GeV/c and p-Be interactions at 
8.9 GeV/c, useful for the understanding of the accelerator neutrino fluxes in the K2K, MiniBooNE and SciBooNE 
experiments, have been published in references [151 116j. 

In this paper we address another of the main motivations of the HARP experiment: the measurement of the yields 
of positive and negative pions relevant for a precise calculation of the atmospheric neutrino fluxes and improved 
modeling of extensive air showers. We present measurements of the double-differential cross-section (P cr'^ / dpdH. 
for positive and negative pion production (in the kinematic range of momentum 0.5 GeV/c < p^r < 8 GeV/c 
and angle 50 mrad < 9-^ < 250 mrad in the laboratory frame) by protons of 12 GeV/c momentum impinging 
on thin cryogenic N2 and O2 targets of 5.5% and 7.5% nuclear interaction length (Ai), respectively. Results for 
the pion production on a thin carbon target in almost the same kinematic region have already been published 
in |17j . Some of those results will be shown again in this paper with a different binning for comparison (see 
Appendix A). These measurements are performed using the forward spectrometer of the HARP detector. Results 
on the measurement of the double-differential tt^ production cross-section in proton-carbon collisions obtained 
with the HARP large-angle spectrometer (100 MeV/c < < 800 MeV/c and 0.35 rad < 9^ < 2.15 rad) are 
presented in a separate article [12]. 

The existing world data for tt^ production on light targets at low beam momentum (< 25 GeV/c) are mainly 
restricted to beryllium targets and with a limited phase space coverage [HI [iHl [20l [2TJ [22j [23] . The work of 
Eichten et al. 22j has become a widely used standard reference data set. In addition to these data, new results 
from the E910 Collaboration have been recently published [24] . 

Carbon is an isoscalar nucleus as nitrogen and oxygen, so the extrapolation to air is the most straightforward. 
Recently the p-C data at 158 GeV/c provided by the NA49 experiment at CERN SPS in a large acceptance 
range have become available f25^ . Relevant data are expected also from the MIPP experiment at Fermilab [26] . 
We would like to mention that the NA61 experiment [27] took first p-C data at 30 GeV/c in the autumn of 2007. 
The foreseen measurements of importance for astroparticle physics are studies of p-C interactions at incoming 
beam momenta 30 GeV/c, 40 GeV/c, 50 GeV/c and tt^-C interactions at 158 GeV/c and 350 GeV/c. 

It is more difficult for experiments to study P-O2 and P-N2 reactions because cryogenic targets are more 
complicated to handle. The results presented in this paper are the first for this type of targets in this energy 
range. 
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1.1 Experimental apparatus 

The HARP experiment [1] [28] makes use of a large-acceptance spectrometer consisting of a forward and large- 
angle detection system. The HARP detector is shown in Fig. [T] A detailed description of the experimental 
apparatus can be found in Ref. [28j . The forward spectrometer - based on five modules of large area drift 
chambers (NDCl-5) [33] and a dipole magnet complemented by a set of detectors for particle identification 
(PID): a time-of- flight wall (TOFW) [3^, a large Cherenkov detector (CHE) and an electromagnetic calorimeter 
(ECAL) - covers polar angles up to 250 mrad. The muon contamination of the beam is measured with a muon 
identifier consisting of thick iron absorbers and scintillation counters. The large-angle spectrometer - based on 
a Time Projection Chamber (TPC) and Resistive Plate Chambers (RPCs) located inside a solenoidal magnet 
- has a large acceptance in the momentum and angular range for the pions relevant to the production of the 
muons in a neutrino factory. For the analysis described here only the forward spectrometer and the beam 
instrumentation are used. 



Drift Chambers 




Figure 1: Schematic layout of the HARP detector. The convention for the coordinate system is shown in the 
lower-right corner. 

The HARP experiment, located in the T9 beam of the CERN PS, took data in 2001 and 2002. The momentum 
definition of the T9 beam is known with a precision of the order of 1% [3lJ. 

The target is placed inside the inner field cage (IFC) of the TPC, in an assembly that can be moved in and 
out of the solenoid magnet. In the cryogenic target setup used for N2 and O2, the gas was liquefied by thermal 
contact with a bath of helium, compressed to 20 bar and then refrigerated to 13 K by adiabatic expansion. 
The refrigerator system was housed inside a vacuum cryostat (typically 2 x 10~^ bar) ending in a vacuum tube 
containing the target. The target arm had 6 cm diameter, 250 fim thick mylar beam entrance and exit windows, 
and the space separating it from the IFC was flushed with nitrogen gas to avoid condensation. The target used 
for the measurements presented here consisted of a 6 cm long, 3 cm of diameter and 125 fim thick mylar cylinder 
and a curved downstream nose, for an actual target volume of about 75 cm^. 

The thickness of the target is equivalent to about 5.5% Xj (4.84 g/cm^) for N2 and 7.5% A/ (6.85 g/cm^) for 
O2. 



The cooling causes a contraction of the target, which induces an uncertainty on its thickness of the order of 1%. 



2 Data Analysis 
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This is taken into account in the uncertainty on the number of target nuclei in section [2.31 

A set of four multi-wire proportional chambers (MWPCs) measures the position and direction of the incom- 
ing beam particles with an accuracy of Rsl mm in position and «0.2 mrad in angle per projection. A beam 
time-of-flight system (BTOF) measures the time difference of particles over a 21.4 m path-length. It is made of 
two identical scintillation hodoscopes, TOFA and TOFB (originally built for the NA52 experiment [32]), which, 
together with a small target-defining trigger counter (TDS, also used for the trigger), provide particle identifi- 
cation at low energies. This provides separation of pions, kaons and protons up to 5 GeV/c and determines the 
initial time at the interaction vertex {to). The timing resolution of the combined BTOF system is about 70 ps. 
A system of two N2-filled Cherenkov detectors (BCA and BCB) is used to tag electrons at low energies and 
pions at higher energies. The electron and pion tagging efficiency is found to be close to 100%. At the beam 
energy used for this analysis the Cherenkov counters select all particles lighter than protons, while the BTOF 
is used to reject ions. A set of trigger detectors completes the beam instrumentation. 

The selection of beam protons is performed using the same criteria as described in [15J . A downstream trigger 
in the forward scintillator trigger plane (FTP) was required to record the event, accepting only tracks with a 
trajectory outside the central hole (60 mm) which allows beam particles to pass. 

The length of the accelerator spill is 400 ms with a typical intensity of 15 000 beam particles per spill. The 
average number of events recorded by the data acquisition ranges from 300 to 350 per spill. 

The absolute normalization of the number of incident protons was performed using 'incident-proton' triggers. 
These are triggers where the same selection on the beam particle was applied but no selection on the interaction 
was performed. The rate of this trigger was down-scaled by a factor 64. 

2 Data Analysis 

2.1 Event and particle selection 

A detailed description of the experimental techniques used for data analysis in the HARP forward spectrometer 
can be found in Ref . [HI [33] . 

With respect to our first paper on pion production in p-Al interactions [15j . a number of improvements to the 
analysis techniques and detector simulation have been made. The most important improvements introduced in 
this analysis compared with the one presented in Ref. [15] are: 

— An increase of the track reconstruction efficiency; 

— Better understanding of the momentum scale and resolution of the detector, based on data, which was 
then used to tune the simulation; 

— New particle identification hit selection algorithms both in the TOFW and in the CHE resulting in much 
reduced background and negligible efficiency losses. In the kinematic range of the current analysis the 
pion identification efficiency is about 98%, while the background from mis-identified protons is well below 
1%; 

— Significant increases in Monte Carlo production have also reduced uncertainties from Monte Carlo statistics 
and allowed studies which have reduced certain systematics. 

Further details of these improved analysis techniques can be found in [Tni[I7]. For the cryogenic targets, dedi- 
cated, high statistics Monte Carlo data were produced using an accurate description of the target geometry. 

At the first stage of the analysis a beam particle type is selected using the beam time of flight system (TOF- 
A, TOF-B) and the Cherenkov counters (BCA, BCB) as described in section 11.11 A value of the pulseheight 
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Table 1: Total number of events and selected pions used in the carbon, nitrogen and oxygen thin target analysis 
at 12 GeV/c, and the number of protons on target as calculated from the pre-scaled trigger count. 



Data set 


C 


Na 


Oa 


Total DAQ events 


1062429 


1375780 


246153 


Acc. beam protons forward 


4375230 


5250240 


840256 


interactions 








n" selected with PID 


8179 


14828 


2817 


7r+ selected with PID 


13530 


23748 


4705 



consistent with the pedestal in both beam Cherenkov detectors rejects electrons, pions and kaons. The beam 
TOF system is used to reject ions, such as deuterons, but at 12 GeV/c is not used to separate protons from 
pions. However, we require time measurements in TOF-A, TOF-B and/or TDS to be present which are needed 
for calculating the arrival time of the beam proton at the target. 

The purity of the proton beam is better than 99%, with the main background formed by kaons estimated to be 
0.5%. This contamination is neglected in the analysis. 

Secondary track selection criteria, described in jl7| , are optimized to ensure the quality of momentum recon- 
struction and a clean time-of-flight measurement while maintaining a high reconstruction efficiency. 

The background induced by interactions of beam particles in the materials outside the target is measured by 
taking data without a target in the target holder ("empty target data"). These data are subject to the same 
event and track selection criteria as the standard data sets. 

To take into account this background the number of particles of the observed type (tt^ , ) in the "empty target 
data" are subtracted bin-by-bin (momentum and angular bins) from the number of particles of the same type. 
The uncertainty induced by this method is discussed in section [2T3l and labeled "empty target subtraction". The 
event statistics is summarised in Table [TJ 



2.2 Cross-section calculation 



The cross-section is calculated as follows 



p'6.,a' 



where 



^;^{Pi-,(^j) is the cross-section in mb/(GeV/c sr) for the particle type a (p, 7r+or 7r~) for each true 
momentum and angle bin (p^, Oj) covered in this analysis; 

N°' (p^, 9'j) is the number of particles of type a in bins of reconstructed momentum p'^ and angle 9'^ in the 
raw data; 

^'p0ap'e'a' is the correction matrix which accounts for efficiency and resolution of the detector; 
jf^jf: Ap^An normalization factors, namely: 

A 



is the number of target nuclei per unit area F^: 
-^pot is the number of incident beam particles on target (particles on target 
Api and AQj arc the bin sizes in momentum and solid angle, respectively 



2 



^' y4 - atomic mass, - Avogadro number, p - target density and t - target thickness 
2) Api = p™'"' - p™'", AQj = 27r(cos{e™'") - cos{e™="')) 



2.2 Cross-section calculation 
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We do not make a correction for the attenuation of the proton beam in the target, so that strictly speaking the 
cross-sections are vahd for a Ai = 5.5% (7.5%) N2 (O2) target. 

The calculation of the correction matrix M^°g ,^^,g,^, is a rather difficult task. Various techniques are described in 
the literature to obtain this matrix. As discussed in Ref. [THj for the p-Al analysis of HARP data at 12.9 GeV/c, 
two complementary analyses have been performed to cross-check internal consistency and possible biases in the 
respective procedures. A comparison of both analyses shows that the results are consistent within the overall 
systematic error jl5j . 

In the first method - called "Atlantic" in [15] - the correction matrix M^°g ,^^,g,^, is decomposed into distinct 
independent contributions, which are computed mostly using the data themselves. The second method - called 
"UFO" in [TS] - is the unfolding method introduced by D'Agostini [33]. It is based on the Bayesian unfolding 
technique. In this case a simultaneous (three dimensional) unfolding of momentum p, angle 9 and particle type 
a is performed. The correction matrix is computed using a Monte Carlo simulation. This method has been used 
in the recent HARP publications [H] [121 [13] and it is also applied in the analysis described here (see [T71 [36] 
for additional information). 

The Monte Carlo simulation of the HARP setup is based on GEANT4 [37] . The detector materials are accurately 
described in this simulation as well as the relevant features of the detector response and the digitization process. 
All relevant physics processes are considered, including multiple scattering, energy loss, absorption and re- 
interactions. The simulation is independent of the beam particle type because it only generates for each event 
exactly one secondary particle of a specific particle type inside the target material and propagates it through 
the complete detector. A small difference (at the few percent level) is observed between the efficiency calculated 
for events simulated with the single-particle Monte Carlo and with a simulation using a multi-particle hadron- 
production model. A similar difference is seen between the single-particle Monte Carlo and the efficiencies 
measured directly from the data. A momentum-dependent correction factor determined using the efficiency 
measured with the data is applied to take this into account. The track reconstruction used in this analysis and 
the simulation are identical to the ones used for the tt"*" production in p-Be collisions [TBI . A detailed description 
of the corrections and their magnitude can be found there. 

The reconstruction efficiency (inside the geometrical acceptance) is larger than 95% above 1.5 GeV/c and drops 
to 80% at 0.5 GeV/c. The requirement of a match with a TOFW hit has an efficiency between 90% and 95% 
independent of momentum. The electron veto rejects about 1% of the pions and protons below 3 GeV/c with 
a remaining background of less than 0.5%. Below Cherenkov threshold the TOFW separates pions and protons 
with negligible background and an efficiency of w98% for pions. Above Cherenkov threshold the efficiency for 
pions is greater than 99% with only 1.5% of the protons mis- identified as a pion. The kaon background in the 
pion spectra is smaller than 1%. 

The absorption and decay of particles is simulated by the Monte Carlo. The generated single particle can re- 
interact and produce background particles by hadronic or electromagnetic processes, thus giving rise to tracks 
in the dipole spectrometer. In such cases also the additional measurements are entered into the migration 
matrix thereby taking into account the combined effect of the generated particle and any secondaries it creates. 
The absorption correction is on average 20%, approximately independent of momentum. Uncertainties in the 
absorption of secondaries in the dipole spectrometer material are taken into account by a variation of 10% of 
this effect in the simulation. The effect of pion decay is treated in the same way as the absorption and is 20% 
at 500 MeV/c and negligible at 3 GeV/c. 

The uncertainty in the production of background due to tertiary particles is larger. The average correction is 
~10% and up to 20% at 1 GeV/c. The correction includes reinteractions in the detector material as well as a 
small component coming from reinteractions in the target. The validity of the generators used in the simulation 
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was checked by an analysis of HARP data with incoming protons, and charged pions on aluminium and carbon 
targets at lower momenta (3 GeV/c and 5 GeV/c). A 30% variation of the secondary production was applied. 
The average empty-target subtraction amounts to «20%. 

Owing to the redundancy of the tracking system downstream of the target the detection efficiency is very 
robust under the usual variations of the detector performance during the long data taking periods. Since the 
momentum is reconstructed without making use of the upstream drift chamber module (which is more sensitive 
in its performance to the beam intensity) the reconstruction efficiency is uniquely determined by the downstream 
system. No variation of the overall efficiency has been observed. The performance of the TOFW and CHE system 
have been monitored to be constant for the data taking periods used in this analysis. The calibration of the 
detectors was performed on a day-by-day basis. 

2.3 Error estimation 

The total statistical error of the corrected data is composed of the statistical error of the raw data and of the 
statistical error of the unfolding procedure, as the unfolding matrix is obtained from the data themselves, thus 
contributing also to the statistical error. The statistical error provided by the unfolding program is equivalent 
to the propagated statistical error of the raw data. In order to calculate the statistical error of the unfolding 
procedure a separate analysis is applied, as described in [17[ I38| . Its conclusion is that the statistical error 
provided by the unfolding procedure has to be multiplied globally by a factor of 2, which is done for the 
analyses described here. This factor is somewhat dependent on the shape of the distributions. For example a 
value 1.7 was found for the analysis reported in Ref. [XT'. 

Different types of sources induce systematic errors for the analysis described here: track yield corrections (~ 5%), 
particle identification (~ 0.1%), momentum and angular reconstruction (~ 0.5%) [f^. The strategy to calculate 
these systematic errors and the different methods used for their evaluation are described in [ITj. An additional 
source of error is due to misidentified secondary kaons, which are not considered in the particle identification 
method used for this analysis and are subtracted on the basis of a Monte Carlo simulation, as in [T7]. No explicit 
correction is made for pions coming from decays of other particles created in the target, as they give a very 
small contribution according to the selection criteria applied in the analysis. 

As a result of these systematic error studies, each error source can be represented by a covariance matrix. The 
sum of these matrices describes the total systematic error, as explained in [17]. 

On average the total integrated systematic error is around 5 — 6%, with a differential bin to bin systematic error 
of the order of 10 — 11%, to be compared with a statistical integrated (bin-to-bin differential) error of ~ 2 — 3% 
(~ 10 — 13%). Systematic and statistical errors are roughly of the same order. 

The overall normalization of the results is calculated relative to the number of incident beam particles accepted 
by the selection. The uncertainty is 2% for incident protons. The contraction of the target with cooling induces 
an additional systematic error of 1% on the N2 and O2 data. 

3 Results 

In Figure [H the measured 7r+ and tt~ spectra in P-N2 and P-O2 interactions at 12 GeV/c are compared to 
an empirical prametrization, developed by Sanford and Wang |39j to describe the production cross-sections of 
mesons in proton-nucleus interactions. The parameters fitted to our p-C data at 12 GeV/c in [T^ have been 

^) The quoted error in parenthesis refers to fractional error of the integrated cross-section in the kinematic range covered by the 
HARP experiment 



3 Results 
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HARP P-N2 -> HARP p-02 -> 77+" 




p [GeV/c] p [GeV/c] 



Figure 2: Measurement of the double-differential production cross-section of positive (open circles) and negative 
(filled circles) pions from 12 GeV/c protons on N2 (left) and O2 (right) as a function of pion momentum, p, 
in bins of pion angle, 0, in the laboratory frame. The curves show the Sanford-Wang parametrization with the 
parameters given in Ref. [T7j (solid line for 7r+ and dashed line for 7r~), computed for the central value of each 
angular bin. In the top right corner of each plot the covered angular range is shown in mrad. 



used and only a constant overall rescaling factor accounting for the target atomic mass has been applied. One 
can observe that the shape and normalization obtained using the carbon data predict quite well the nitrogen 
and oxygen data. This point will be made more clear when the N2/C and O2/C ratios are taken. The shapes 
of the momentum spectra are similar for secondary 7r+ and tt" , as well as for different data sets, where only a 
different normalisation factor can be noticed because of the different nuclear masses of the target nuclei. The 
conclusions drawn in |17j appear to be confirmed for the data sets presented here: the parametrization provides 
an approximate description of the main features, but is not able to describe the data well in some regions of 
kinematic space, particularly at high momenta and at large angles. 

The central values and square-root of the diagonal elements of the covariance matrix are listed in Tables [2] 
and [31 

The kinematic range of the measurements covers the momentum region from 0.5 GeV/c to 8.0 GeV/c (subdivided 
into 10 intervals) and the angular range from 0.05 rad to 0.25 rad (subdivided into 4 intervals). The error 
bars correspond to the combined statistical and systematic errors as described in section 12.31 The overall 
normalization error of 2% for the normalization of incident protons and of 1% for the target size variation are 
not shown. 

The pion production ratios N2/C and O2/C are presented in Figs. [3]l4] and are compared to GEANT4 Monte 
Carlo predictions. As noted before, the difference between the target materials is justified by an overall nor- 
malisation factor taking into account the different nuclear masses of the target materials. The various models 
(see [40] for details) do predict the ratio of cross-sections accurately, with very little spread between them. This 
conclusion is different when the absolute predictions models are compared with the measured cross-sections as 
shown in Ref. [T7] . 
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Table 2: HARP results for the double-differential tt"*" and tt" production cross-section in the laboratory system, 
(Pa'^ /{dpdCl), for p-N2 interactions at 12 GeV/c. Each row refers to a different (Pmin ^ P < Pmax, ^min ^ ^ < 
^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



nil 11 
(rad) 


Omax 

(ra.,1) 


Pmin 

(Ge\-/c) 


Pmax 
(GcV/c) 


(Pa''^ /{(IpdO.) 
(ml)/'(CT('\'/'(: sr)) 


(ml)/(C, 


/(rfpdO) 
o\-/c sr)) 


0.05 


0.10 


0.50 


1.00 


240.5 


± 


39.9 


200.7 


± 


33.7 






1.00 


1.50 


318.1 


± 


29.0 


256.9 


± 


27.9 






1.50 


2.00 


288.2 


± 


28.1 


195.9 


± 


23.2 






2.00 


2.50 


310.5 


± 


26.3 


165.3 


± 


18.0 






2.50 


3.00 


280.6 


± 


23.9 


127.4 


± 


14.8 






3.00 


3.50 


183.0 


± 


18.9 


71.2 


± 


10.1 






3.50 


4.00 


152.3 


± 


15.6 


64.7 


± 


8.8 






4.00 


5.00 


98.3 


± 


9.5 


37.5 


± 


5.6 






5.00 


6.50 


46.4 


± 


5.8 


14.3 


± 


3.0 






6.50 


8.00 


14.4 


± 


2.7 


2.4 


± 


0.8 


0.10 


0.15 


0.50 


1.00 


239.5 


± 


35.5 


304.3 


± 


42.8 






1.00 


1.50 


336.9 


± 


34.7 


217.3 


± 


24.5 






1.50 


2.00 


250.0 


± 


27.2 


180.8 


± 


20.9 






2.00 


2.50 


229.3 


± 


25.6 


83.9 


± 


15.4 






2.50 


3.00 


145.3 


± 


17.5 


78.7 


± 


11.9 






3.00 


3.50 


101.8 


± 


13.4 


63.0 


± 


8.6 






3.50 


4.00 


61.9 


± 


9.6 


39.3 


± 


6.5 






4.00 


5.00 


33.9 


± 


5.1 


14.3 


± 


3.4 






5.00 


6.50 


11.3 


± 


2.2 


3.7 


± 


1.0 






6.50 


8.00 


2.0 


± 


0.6 


0.5 


± 


0.2 


0.15 


0.20 


0.50 


1.00 


337.6 


± 


45.5 


269.8 


± 


40.0 






1.00 


1.50 


244.4 


± 


27.8 


151.1 


± 


21.8 






1.50 


2.00 


148.6 


± 


19.3 


122.4 


± 


18.7 






2.00 


2.50 


107.4 


± 


15.8 


83.7 


± 


14.7 






2.50 


3.00 


75.8 


± 


13.6 


42.3 


± 


10.5 






3.00 


3.50 


36.5 


± 


7.7 


17.5 


± 


6.7 






3.50 


4.00 


27.7 


± 


5.2 


20.3 


± 


5.2 






4.00 


5.00 


17.6 


± 


3.1 


6.0 


± 


2.1 






5.00 


6.50 


4.9 


± 


1.5 


0.4 


± 


0.3 






6.50 


8.00 


0.6 


± 


0.2 








0.20 


0.25 


0.50 


1.00 


212.6 


± 


36.0 


145.5 


± 


28.0 






1.00 


1.50 


84.1 


± 


20.3 


193.4 


± 


34.2 






1.50 


2.00 


95.4 


± 


22.3 


123.2 


± 


25.2 






2.00 


2.50 


63.2 


± 


14.2 


84.4 


± 


18.5 






2.50 


3.00 


47.6 


± 


10.9 


49.8 


± 


11.8 






3.00 


3.50 


21.7 


± 


5.9 


15.5 


± 


5.2 






3.50 


4.00 


9.2 


± 


3.6 


6.2 


± 


4.2 






4.00 


5.00 


7.7 


± 


2.6 


1.5 


± 


2.2 






5.00 


6.50 


1.9 


± 


1.3 












6.50 


8.00 
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Table 3: HARP results for the double-differential tt"*" and tt" production cross-section in the laboratory system, 
(Pa'^ /{dpdCl), for P-O2 interactions at 12 GeV/c. Each row refers to a different (Pmin < Pmax, ^min ^ ^ < 
^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



nil 11 


Omax 


Pmin 


Pmax 




^/(rfjxin) 




/(rfpdO) 


(rad) 


(ra.,1) 


(Ge\-/c) 


(GcV/c) 


(ml)/( 


Gc\' 




(ml)/(C, 


v\'/ 

1 


c sr)) 


0.05 


0.10 


0.50 


1.00 


290.3 


± 


63.9 


290.6 


± 


63.3 






1.00 


1.50 


417.5 


± 


56.5 


307.7 


± 


49.4 






1.50 


2.00 


435.2 


± 


53.9 


236.9 


± 


41.2 






2.00 


2.50 


371.4 


± 


47.3 


166.1 


± 


31.4 






2.50 


3.00 


313.9 


± 


42.0 


140.5 


± 


27.3 






3.00 


3.50 


266.2 


± 


36.9 


117.5 


± 


23.3 






3.50 


4.00 


175.3 


± 


27.1 


71.1 


± 


16.7 






4.00 


5.00 


127.4 


± 


18.1 


41.7 


± 


9.9 






5.00 


6.50 


63.0 


± 


10.3 


20.8 


± 


6.1 






6.50 


8.00 


25.0 


± 


6.8 


7.3 


± 


3.3 


0.10 


0.15 


0.50 


1.00 


327.5 


± 


62.5 


462.9 


± 


85.0 






1.00 


1.50 


492.2 


± 


65.7 


330.2 


± 


50.0 






1.50 


2.00 


317.9 


± 


46.3 


235.9 




40.6 






2.00 


2.50 


332.7 


± 


48.5 


107.5 




26.4 






2.50 


3.00 


224.1 


± 


36.6 


138.4 


± 


28.8 






3.00 


3.50 


113.0 


± 


21.9 


72.6 


± 


18.5 






3.50 


4.00 


97.8 


± 


21.7 


57.4 


± 


15.3 






4.00 


5.00 


52.2 


± 


10.9 


24.1 


± 


7.6 






5.00 


6.50 


15.1 


± 


5.4 


3.7 


± 


2.1 






6.50 


8.00 


2.2 


± 


1.1 


1.4 


± 


1.2 


0.15 


0.20 


0.50 


1.00 


337.5 


± 


69.6 


244.2 


± 


56.7 






1.00 


1.50 


313.5 


± 


52.3 


198.7 


± 


42.9 






1.50 


2.00 


225.6 


± 


40.9 


113.9 




30.5 






2.00 


2.50 


144.9 


± 


30.4 


98.7 




28.0 






2.50 


3.00 


78.8 


± 


22.3 


45.5 




19.3 






3.00 


3.50 


55.2 


± 


16.9 


15.2 


± 


10.0 






3.50 


4.00 


29.3 


± 


11.4 


27.5 


± 


12.7 






4.00 


5.00 


19.8 


± 


8.0 


11.0 


± 


6.7 






5.00 


6.50 


6.0 


± 


4.5 












6.50 


8.00 














0.20 


0.25 


0.50 


1.00 


235.0 


± 


59.7 


186.7 


± 


50.2 






1.00 


1.50 


137.0 


± 


43.4 


200.3 


± 


58.2 






1.50 


2.00 


126.2 


± 


42.0 


117.5 


± 


38.2 






2.00 


2.50 


70.2 


± 


31.5 


115.8 


± 


39.6 






2.50 


3.00 


37.3 


± 


16.1 


60.8 


± 


25.7 






3.00 


3.50 


36.8 


± 


16.7 


22.1 


± 


15.0 






3.50 


4.00 


18.0 


± 


13.2 


15.1 


± 


15.3 






4.00 


5.00 


11.0 


± 


9.2 


4.5 


± 


5.5 






5.00 


6.50 


2.8 


± 


4.8 












6.50 


8.00 















4 Summary and conclusions 

The results reported in this article may contribute to the precise calculations of atmospheric neutrino fluxes and 
to the improvement of our understanding of extensive air showers simulations and hadronic interactions at low 
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Figure 3: p-N2 to p-C production ratio for 7r+(left panel) and tt^ (right panel) at 12 GeV/c, compared with 
GEANT4 simulation predictions using different models. In the top right corner of each plot the covered angular 
range is shown in mrad. Only statistical errors are used, since most systematic ones cancel. 



energies. A detailed study of the role of hadronic interactions for production of muons in extensive air showers, 
which are one of the main ingredients to infer the mass and the energy of the primary cosmic ray particle, is 
shown in reference [3S] . 

In this paper we presented measurements of the double-differential production cross-section of positive and 
negative pions in the collisions of 12 GeV/c protons with thin nitrogen and oxygen targets. The data were 
reported in bins of pion momentum and angle in the kinematic range 0.5 GeV/c < Ptt < 8 GeV/c and 0.05 rad 
< 6„ < 0.25 rad in the laboratory frame. A detailed error analysis has been performed yielding total bin-to-bin 
differential errors (statistical and systematic) of about 15%, an overall normalization error of 2% and additional 
1% for the target size variation. We should stress that the HARP data are the first measurements with cryogenic 
targets in this kinematic region with good precision. 

Simulations predict that collisions of protons with a carbon target are very similar to proton interactions with 
air (see e.g. [33]). That explains why these datasets can be used for tuning models needed in astroparticle physics 
simulations. Our measurements on p-N2 and P-O2 confirm these predictions. 
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A Cross-section data 

The tabulated cross-section data for p-C interactions at 12 GeV/c, already published in [TT], are reported again 
here with a different binning for comparison. 

Table 4: HARP results for the double-differential 7r+ and tt" production cross-section in the laboratory system, 
(Pa'^ /{dpd^), for p-C interactions at 12 GeV/c. Each row refers to a different (Pmin ^ P < Pmax, ^min < ^ < 
^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 
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